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Abstract 

In the minimal super symmetric standard model, large left-right mixing of stau's 
is sometimes intriguing from phenomenological viewpoints. However, too large 
left-right mixing is not acceptable since the electroweak-breaking vacuum becomes 
metastable. In this paper the vacuum transition rate is evaluated by using semi- 
classical techniques, and constraints on parameters of the model are shown. In the 
calculation the bounce solution is derived in the multifield space. These constraints 
are also applied to the case of the low-energy minimal gauge mediation model. 



1 Introduction 



The minimal supersymmetric standard model (MSSM) has two Higgs doublets, and the 
ratio of their vacuum expectation values, tan (3(= (H u ) / (H^)), is an important parameter 
when the phenomenology is discussed. It is expected to be between ~ 2 and ~ 60 if the 
Yukawa coupling constants for the third generation fermions are perturbative below the 
GUT scale. 

Large tan /3 is favored from several phenomenological viewpoints. First, the light Higgs 
boson mass depends on tan with a high mass at larger tan /3. The muon (g — 2) anomaly 
also favors large tan/3 since the SUSY correction is proportional to tan/3 pQ. Second, the 
Yukawa unification of the third generation fermions in the SO (10) SUSY GUT requires 
large tan (3 [2J. 

Third, the SUSY CP problem is automatically solved when A and _B M parameters 
vanish at high energy scale. The A and parameters are for trilinear and Higgs bilinear 
soft-SUSY breaking couplings, respectively. In the case, large tan /3 is predicted, because 
tan/3 is proportional to inverse of the parameter, while the B^ parameter is only 
radiatively generated. It is pointed out in Refs. [3J H] that in the low-energy minimal 
gauge mediation (MGM) model, in which the messenger scale is around lO*- 5-6 -' GeV, 
tan/3 is predicted to be (50 — 60). The radiative correction induces tiny B^ parameter 
due to the low messenger scale. 

Forth, in some gauge mediation models, stau is the next-lightest SUSY particle (NLSP). 
It is so long- loved that the big bang nucleosynthesis (BBN) may be destroyed [5] . When 
the left-right mixing term in the stau mass matrix, which is proportional to tan /3, is large, 
the annihilation of stau pair is enhanced due to the Higgs s-channel exchange so that the 
primordial abundance of stau NLSP is reduced [6]. 

While large tan/3 sometimes intriguing as mentioned above, it is known that large 
tan/3 solutions may suffer from vacuum instability |4J. When the left-right mixing term 
for stau's is increased, electric charge-breaking minimum in the scalar potential appears 
and it becomes deeper than than the "ordinary" electroweak-breaking minimum. The 
lifetime of the electroweak-breaking vacuum is required to be longer than the age of the 
universe. 

In this paper, we derive upperbound on the left-right mixing term for stau's as a 
function of the left- and right-handed slepton masses, by imposing that the electroweak- 
breaking vacuum has longer lifetime than the age of the universe. The corresponding 
lowerbound on the stau mass is also shown as a function of the left- and right-handed 
slepton masses. It is found that the constraints are insensitive to the parameters in the 
Higgs potential, including tan/3 itself. When the LHC experiment discovers stau's or 
gives lowerbound on the mass, the stability of the electroweak-breaking vacuum would be 
useful to derive the parameters in the stau mass matrix. 

The quantum transition rate of the metastable vacuum is estimated by semiclassical 
technique [TJ. We evaluate numerically bounce configuration in three- fields space (left- 
and right-handed stau's and Higgs boson). The quantum transition rate is evaluated in the 
previous works, though the bounce configuration is approximated to be one-dimensional 
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[U E]. Thus, our result is more accurate than those. 

This paper is organized as follows. In next section, the quantum transition rate of the 
electroweak-breaking vacuum to the true one is evaluated, and stability of the electroweak- 
breaking vacuum is discussed. The constraints on parameters for stau mass matrix are 
derived there. In Section 3 the constraints on the low-energy MGM model are shown. 
Section 4 is devoted to conclusions and discussion. 



2 Stability of vacuum 

Stability of the vacuum in the MSSM in the case that tan 8 is large have been studied in 
Ref. jl]. In this section we follow their argument and give a more reliable evaluation for 
vacuum stability. 

In gauge mediation models and the constrained MSSM which demands unification 
of soft scalar squared masses at high scale, sleptons are lighter than squarks because in 
the former cases slepton masses are proportional to small squared gauge couplings g Y of 
U(l)y and g\ of SU(2) L and in the latter case squarks receive large gluino contribution 
through renormalization group (RG) running. If tan 8 is large the lightest sfermion is 
stau for large left-right mixing, which is proportional to the Yukawa coupling of the tau 
lepton, 

m T tan 8 , , 

y v d 100 V ; 

Too large tan 8 makes stau tachyonic and classical stability of the electroweak-breaking 
vacuum is lost. Although the electroweak-breaking vacuum is classically stable, it may be 
unstable by quantum tunneling effect. If the electroweak-breaking vacuum is not a global 
minimum but local one in the potential, it would collapse into the global one eventually. 
Such a metastable vacuum is viable only when its lifetime is longer than the age of the 
universe. 

To see the situation described above in detail we write the scalar potential for neutral 
component of up- type Higgs H u , left-handed stau L and right-handed stau f R as follows, 

V = {m 2 Hu + fi 2 )\H u \ 2 + m 2 L \L\ 2 + m 2 fR \f R \ 2 - {y Tf iH*Lf R + h.c.) + y 2 T \L 2 f R \ 2 

+ j(\L\ 2 + l^l 2 ) 2 + 9 f(\L\ 2 ~ 2\tr\ 2 - \H u \ 2 f + g -^^S H \H u \\ ^ 

where /i is supersymmetric Higgs mass and m 2 's are soft squared masses of each scalar. 
For simplicity, /i is assumed to be real. Terms including down-type Higgs are ignored 
since the vacuum expectation value (VEV) is small for large tan 8. The scalar potential 
includes an RG-improved term, which is the last term in Eq. ([2]) and reflects quartic Higgs 
interaction induced by a loop diagram of the top quark: 

*H = - 2 ^— 2 log^. (3) 

7T 2 9y + 92 m t 
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Here, y t is the Yukawa coupling of the top quark and m t and m ( - are the mass of the top 
and stop, respectively. Typical value of 5# is ~ 1. In addition, loop diagrams of the stop 
also induces the quartic interaction and may make sizable contribution in some cases. It 
is proportional to Af, where A t is an A parameter for stop interaction. 

The H u mass term, m 2 H + /i 2 , is negative so that H u gets a VEV v u . Expanding 
around this electroweak-breaking vacuum, we have 



V= m W( m l + 9 -t^tvl)L> + (rnl R 
- 2y T fj,(j)LT R + 2 



















R+ v u 


<f> 3 







(4) 



where H u = v u + and m 2 = gy ^ ga (1 + fe)f 2 = (1 + 5h) sin 2 /3m|. Ellipsis stands for 
quartic terms. We show only terms with the real parts of scalar bosons in above scalar 
potential. The first line of Eq. (TjJ gives a squared mass matrix of stau's, 

M 2 ~ = ( m L + ^~ S w) m Z WtVu \ , 5 n 

V WtVu m 2 ~ R + s 2 w m 2 z J ' 

Here s^y = gy/igy + g^) an d we take a large tan/3 limit. Classical stability of the 
electroweak-breaking vacuum is equivalent to the positivity of the smaller eigenvalue of 
the matrix Aif. 

Even if the classical stability condition is satisfied, the first term in the second line 
of Eq. could generate a global minimum where (0), (L), and (f R ) ^ and make the 
electroweak-breaking vacuum metastable. 

Quantum transition rate of the metastable vacuum is estimated by semiclassical tech- 
nique [7j. In this technique the imaginary part of energy of the false vacuum, which 
is proportional to the transition rate is evaluated using path integral method in Eu- 
clidean spacetime. The path integral is dominated by so-called bounce configuration, 
(Pi = (fi(t,x), where (fx — (p, (p 2 — L and <£> 3 = f R in our case. It is a stationary point of 
the action and satisfies boundary conditions \im t ^,± OQ {p i (t,x) = ip{, where <p{ are values 
of the fields at false vacuum. It is known that we can take an 0(4) symmetric solution 
[8]. The 0(4) symmetric Euclidean action is 



(6) 



where r is a radial coordinate in four- dimensional spacetime. The equation of motion and 
boundary conditions are 

d?<Pi 6dipi dV 
Km <Pi(r) = V {, ^(0) = 0. (8) 
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Figure 1: Bounce configuration <£>j(r) a s a function of the radial coordinate r. Parameters 
of potential and B are shown in figure. 



Finally the vacuum transition rate per unit volume is evaluated as follows, 

T/V = Ae- B . (9) 

The prefactor A is the fourth power of the typical scale in the potential. Its precise value 
is hard to calculate but the transition rate is not sensitive to it. On the other hand, B 
has an importance in the evaluation and 

B = S E [<p(r)]-S E [(p f ). (10) 

If we demand that T/V is much smaller than the fourth power of the present Hubble 
expansion rate Hq = 1.5 x 10 -42 GeV and assume A is (100 GeV) 4 , we have constraint 
on metastability that B > 400. 

For numerical calculation of bounce solutions we used the method of Ref. [9]. Search 
for bounce solution with single field is easily achieved by overshooting/undershooting 
method in which one scans initial values <f{0) around the global minimum. This scan is 
difficult in the case of multifield. The authors of Ref. j9] found that in one-dimensional 
space the initial value <p(0) is obtained by using a modified potential. Once we have a 
bounce solution in one-dimensional space, gradually increasing the spacetime dimension 
we obtain the four-dimensional solution. Actual calculation is based on a discretized 
equation of motion. 

This method is unsuitable yet to the case that the global minimum is much deeper than 
local one. Because in that case the modified potential is quite different from original one 
so that it is difficult to obtain the one- dimensional solution. For such parameter regions, 
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starting from the solution obtained in region where two minima are nearly degenerate, 
we gradually change the parameters of the potential and reach the solution. 

Fig. [I] is a bounce solution obtained by this method. A trajectory is slightly different 
from a straight line assumed in Ref. [I]. When ~ 0.2v u , L and approach almost 0. 

For the potential of Eq. (0J relevant parameters to vacuum stability are m|, rn~ R , //, 
tan/3 and Sh- Since the coefficient of quadratic term Ltr is proportional to //tan/3 and 
that of cubic term 4>Ltr is nearly proportional to //tan/3, the bounce solution is sensitive 
to //tan/3, but not tan/3 itself. To confirm this, we plot contours of B changing //tan/3 
and tan/3 while other parameters are fixed (Fig. [2] (top)). Irrelevance of tan /3 is clearly 
seen in the figure. Weak dependence on tan /3 is generated by the quartic term xfcl?f\ in 
Eq. (|2J). This term lifts the global minimum and increases B a little. 

Next, we also plot contours changing //tan/3 and Sh in Fig. [2] (bottom). When A t 
is negligible, Sh mainly depends on the stop mass. If we take in the range 600 GeV 
< mi < 2 TeV, Sh has the value from 0.65 to 1.3. Larger value of Sh stabilizes the 
electroweak-breaking vacuum and slightly increases B. Thus, it is found that the transi- 
tion rate is sensitive to only m|, fn\ R , and //tan/3 among parameters in the potential of 
Eq. OH). 

In Fig. |3l we show upperbound on // tan /3 that satisfies the metastability condition 
B > 400. In numerical calculation we set Sh — 1 and tan/3 = 70. It is found from 
Fig. [2] (top) that dependence on tan/3 changes the result at most 5%. We also evaluate 
that the dependence on Sh is less than 1% in the range that 0.65 < Sh < 1.3. We also 
show upperbound on // tan (3 that satisfies a rigorous stability condition (the electroweak- 
breaking vacuum is global minimal in the potential) in Fig. HI The metastability condition 
gives about 50% looser constraint on // tan (3 than that the rigorous stability condition. 

Upperbound on // tan (3 gives lowerbound on stau mass. We plot stau mass and mixing 
angle when // tan /3 is taken to be its maximum value allowed by metastability in Fig. 
Since LEP2 experiments give lowerbounds on stau mass as m fl > 80-90 GeV [TU], in 
these region in Fig. [5] the metastability condition gives severer constraint on //tan/3 than 
the experimental bounds. 

We fit this result under the metastability condition to a function as follows, 

//tan/3 < 76.9y/mim,T R + 38.7(m z + m fR ) - 1.04 x 10 4 GeV. (11) 

The difference of this fit is less than 1% in the region where > 80 GeV. This would 
be useful when considering large left-right mixing of stau's. 

3 Application 

As an application of the previous section, we consider the low-energy MGM model. The 
MGM model is the gauge mediation model with a constraint that the B^ term vanishes 
at the messenger scale. This model is a solution for the SUSY CP problem in Refs. [21 H] , 
as mentioned in Introduction. 
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In this model tan /3 is an output of the model rather than input. At the tree level, 
and tan j3 are related as follows, 

2jE?At -sin 2/3 (12) 



m 2 Hu + m 2 Hd + 2/i 2 ~ tan/3 

Since is generated only through RG running, large tan f3 is predicted in the low-energy 
MGM model. 

There are two classes of contributions to B^. The first one arises from one- loop 
diagrams of gaugino. The other class is generated by A terms. Since A terms also vanish 
at the messenger scale, this contribution arises at two loop level. The largest one among 
A term contributions is proportional to squared Yukawa coupling of the top (or bottom) 
quark and squared SU(2>)c gauge coupling. Because of these large couplings, the A term 
contributions are comparable with the gaugino one despite the extra loop. These two 
contributions have opposite signs, and B^ is smaller than naive expectation in the low- 
energy MGM model. 

We assume that messengers are 5 + 5* in SU(5) and write couplings to SUSY-breaking 
sector field S as follows, 

W = X T Sqq + X D Sll, (13) 

where, q and / are SU(3)c triplet and SU(2)l doublet in messenger, respectively. A's are 
coupling constants and S has VEVs (S) = M + 2 F$. When A's are equal at the GUT 
scale, they differ at messenger scale due to RG running, which takes At larger than Xd 
at low energy. The gaugino masses are expressed as follows, 

M, = (14) 
where oti = g 2 /^ and the loop function g is 

g{x) = — [(1 + x) log(l + x) + (1 - x) log(l - x)]. (15) 

x 2) 3 are defined by x 2 = F s /XdM 2 , x 3 = F S /X T M 2 and g(x\) = \g{x^) + |g(x 2 )- The 
difference in A's affects the gaugino masses only through the loop function and its effect 
becomes large in the low-energy MGM model where F$ ~ M 2 . 

In Fig. E] we show the prediction of the low-energy MGM model and constraints by 
(meta-)stability of the vacuum. We take Xd = 1.3Xt (top), Xt = Xp (middle), At = 1.3Ao 
(bottom) and max^,^) = 0.95. These results are obtained using SuSpect |TT] with 
modification for splitting in couplings of messenger. In the top and bottom figures the ratio 
of loop functions is #(0.95) : #(0.95/1.3) = 1.14 : 1. In gauge mediation models current 
experimental lowerbound on stau mass is 87.4 GeV [10]. In the case of At = 1.3Xd, tan/3 
is larger than 60 and the metastability condition gives severer constraint on lower limit 
of M2 than experimental lowerbound on stau mass gives. This is explained as follows. In 
Fig. |6]the magnitude of the A-term contributions to B^ is larger than that of the gaugino 
contributions. If At > Ad, the gluino mass is smaller than the value that GUT relation 
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is valid. Since leading parts in A-term contribution is proportional to gluino mass, the 
decrease in gluino mass enhances the cancellation between two classes of contributions. 
As a result, _B M becomes smaller at the electroweak scale and larger tan/3 is predicted. 
Our result is consistent with those in Ref. [I]. 

4 Conclusions and discussion 

In this paper we have studied vacuum stability constraints on left-right mixing of stau's. 
/xtan/3 is bound by quantum transition to the charge-breaking vacuum. We have shown 
lowerbound on stau mass from the constraint and its corresponding mixing angle. Our 
estimation to vacuum transition rate is based on semiclassical technique performed in the 
three-fields space, which give more precise results than previous works. 

We have also considered these constraints in the low-energy MGM model. The vacuum 
stability conditions constrain parameters of model more severely than experimental limit 
of stau mass in the case that messenger triplet is heavier than doublet. Such split of 
couplings in messenger is expected from RG running. 

One of issues that we have not handled in this paper is the vacuum transition in the hot 
universe. The condition that thermal transition to charge-breaking vacuum is suppressed 
enough would also restrict /xtan/3. In addition, electroweak symmetric vacuum should not 
decay to the charge-breaking vacuum at finite temperature. To trace electroweak-breaking 
process properly is beyond the scope of this paper and will be studied elsewhere. 
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Figure 2: Solid lines are contours of B = 100, 400 and 1000 in tan/3-/xtan/3 plane (top) 
and Sh~ Aitan/3 plane (bottom). The electroweak-breaking vacuum becomes metastable 
above broken lines. 
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Figure 4: Contour plots for upperbound on //tan/3 that satisfies condition that the 
electroweak-breaking vacuum is global minimum. 



10 



100 



100 200 300 400 

7Uf R [GeV] 



500 



600 



Figure 5: Contour plots for stau mass (solid lines) and stau mixing angle Of (broken 
lines) in m^-m^ plane when /i tan f3 is taken to be upperbound shown in Fig. [3j 



11 



Ap/ At — 1.3 



70 
60 
50 
40 
30 
20 



unstable 



B > 400 

MGM prediction 



metastatic 



/ m f . > 87.4 GcV 



100 150 200 250 300 350 400 450 500 
M 2 [GcV] 



At — Ar 




100 150 200 250 300 350 400 450 500 

M 2 [GeV] 

\t i = 1-3 




100 150 200 250 300 350 400 450 500 

M 2 [GeV] 

Figure 6: (Meta-)stability constraints in the low-energy MGM model. Horizontal axis 
corresponds to gaugino mass of SU(2)l, and vertical axis corresponds to tan/3. Heavy 
lines indicate predicted value of tan (3 in MGM. Electroweak-breaking vacuum is unstable 
above broken lines and metastable above dotted lines. Below solid lines metastability 
condition B > 400 is satisfied. In the right-hand side of dot-dashed lines stau mass is 
larger than experimental lowerbound. We take \d = 1.3 Ay (top), At = A^ (middle), 
At = 1.3Ad (bottom) and max(z 2 ,i3) = 0.95 in all cases. 
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